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Abstract. A method of G. Wilson for generating commutative algebras of 
^_ ' ordinary differential operators is extended to higher dimensions. Our construc- 

, tion, based on the theory of ©-modules, leads to a new class of examples of 

commutative rings of partial differential operators with rational spectral vari- 
eties. As an application, we briefly discuss their link to the bispectral problem 
and to the theory of lacunas. 

o' 
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CN ■ 1. Introduction 

1.1. The purpose of this paper is to discuss a certain algebro-geometric construc- 
tion of commutative rings of partial differential operators with singular spectral 

ON ■ 

varieties. The principal idea of the method to be used goes back to the work of 
Krichever (see, e.g., for a review), though our effort here has been inspired by 
f~>) | more recent developments. These are the study of algebraic Schrodinger operators 

initiated by O. Chalykh and A. Veselov |]| and the work of G. Wilson @ || 
C*"** ' on "bispectral" algebras of ordinary differential operators. 

ON i Our attempt is to merge the ideas of the above mentioned authors within a more 

general setting and to explore some new instructive examples. 



1.2. In [|jjj-|16| a certain geometric axiomatics for the Baker-Akhiezer (BA) func- 
tion has been developed in affine spaces of arbitrary dimension (n > 1). The con- 
O • struction involves a finite set A of linear (homogeneous) hyperplanes in C™ with 

prescribed (integer) multiplicities. These data are postulated to define a "pole di- 
visor" of an associated BA function. When it exists, such a function is necessarily 
unique and solves a (formal) spectral problem for a certain commutative ring TZa 
of partial differential operators containing a second order (Schrodinger) operator. 
Algebraically, the ring TZa turns out to be supercomplete in the sense that the min- 
imal number of generators of TZa (as a C-algebra) exceeds its Krull dimension, i.e. 
the dimension of the spectrum Spec7?._4, while geometrically the latter proves to be 
a singular variety. Within the "rational" version of this axiomatics, the differential 
operators in TZa have constant principal symbols and rational lower order coeffi- 
cients with singularities located on the hyperplanes of A. These operators possess a 
number of interesting analytic properties and applications among which one should 
mention the link to Huygens' principle and the theory of lacunas (see || ||, [?]. ||). 

The principal problem within the Veselov- Chalykh approach is to determine all 
possible hyperplane arrangements for which a nontrivial BA function does actually 
exist. Up to now, this problem has been completely settled in dimension n — 2 (see 
§ |), while in higher dimensions (n > 3) only partial results and examples are 



available (see |33|, |3J, |35| ) 



1.3. One may pose a more general question (cf. pS| , Section 4): What algebraic 
(or analytic) sets in C™ are admissible as a singularity locus for the BA function 

l 
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associated to a (non-trivial) supercomplete commutative ring of partial differential 
operators? 

In this paper we give a partial answer to this question in the "rational" (al- 
gebraic) situation. Let X be any open quasi-affme algebraic variety in C™, i.e. 
X = C n \r- x (0) with some t G C[sei, . . . ,x n ], and let O(X) and V(X) stand for 
the rings of regular functions and differential operators on X respectively. It is 



a consequence of Theorem 4.2 below that there exists a nontrivial supercomplete 
commutative ring TZx C D(X) of differential operators on X whose common eigen- 
function ip = t/>(x,£) is unique (up to a normalization) and quasi-regular, i.e. has 
the form ip{x,£) := P(x, £)e( x & with P e 0{T*X) = 0(X)®C[€\. 

This result implies that the BA function (regarded as a function of x only) can 
admit singularities along arbitrary algebraic hypersurfaces in C" . In Section |I] we 
will give a simple and constructive proof of this observation. A number of explicit 
examples will be treated in Section ||. 

It is worth mentioning that the ring TZx does not necessarily contain a sec- 
ond order operator and, hence, our examples go beyond the framework of the 
Veselov-Chalykh axiomatics; however, as we will see in Section @, they have equally 
interesting analytic applications. 

Unlike the BA functions that arise from the Veselov-Chalykh construction, a 
generic ^-function associated with a ring TZx is not self-dual. Nevertheless, we 
demonstrate the existence of a complete ring of partial differential operators in the 
spectral parameter having ip as a common eigenfunction. In this way, we get an 
example of a nontrivial bispectral involution in the case of several variables. 

1.4. For ordinary differential operators, a complete answer to the question stated 
above is a simple consequence of a method of G. Wilson |37j]. His technique was 
reformulated later in terms of Darboux transformations in [jj || and ||2l|, The 
method used in the present work can be regarded as a further generalization of 
this approach with a view to address the problem for partial differential operators 
(n > 1). A natural framework for such a generalization is provided by the theory 
of V-modules. Our intention is to demonstrate that the language of X>-modules is 
quite relevant and illuminating in this context. It brings some new insight even in 
the case n = 1 (see Section ^). In this paper we focus mostly on the construction 
of families of commuting differential operators rather than on the study of their 
properties. In particular, we only briefly discuss the bispectral problem which was 
one of our original motivations. We believe, however, that the understanding of 
this phenomenon within the theory of X>-modules is a promising problem which 
deserves a separate and more detailed investigation. 

2. The theory of ©-modules: basic concepts 

We start with a brief summary of some results and definitions from the theory 
of ©-modules. Our consideration will be restricted to the purely global (algebraic) 
case. For a more general and comprehensive treatment of the subject, the reader 



is referred to the literature |12] [L0|, |27[ [28| (see also |29|, 30 ) 



2.1. Regular Differential Operators. Let A be a non-singular irreducible affinc 
algebraic variety over C. We may (and will) regard X as being imbedded in the 
standard n-dimensional afnne space C", i.e. we identify X with a set of common 
zeros of a (prime) ideal X(X) in the polynomial ring C[x]. Let 0(X) be an affine 
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algebra associated with X (the coordinate ring of regular functions on X). By our 
identification, we have O(X) = C[x]/T(X). Then, O(X) is an integral domain. We 
write T{X) for the quotient field of 0{X). 

Let $f{X) be the set of all C-linear derivations on F{X) and let 0o(X) C Qf(X) 
be its proper subset that preserves the subring O(X) C ^F(X), i.e. qq(X) := {d G 
qAX) I 8 G End c (G(X))}. 

Definition 2.1. The ring of regular (algebraic) differential operators on X is the 
ring of C-linear endomorphisms of 0(X) generated by the derivation operators from 
0O (X) and the multiplication operators from 0{X) itself. This ring is denoted by 
V(X). 

The following examples will be of basic concern in the present work. 

Example 1. If X = C" then V(X) is the Weyl algebra A n (C) = C(x, J^}, i.e. the 
algebra of differential operators with polynomial coefficients. 

Example 2. Let X be any smooth algebraic variety. For r G O(X), 1/r is defined 
as a regular function on the complement X T := X\t^ 1 (Q) of the zero set of r in 
X. We write 0(X) T for the algebra generated by 0{X) and 1/r and call it the 
localization of the ring 0{X) at r. Similarly, the ring V(X T ) of regular differential 
operators on X T can be identified with a localization of T>(X) at r, i.e. 

(1) 2?(X T ) = P(X) r := 0(X) T ® 0{x) V(X). 

In particular, if X = C™ then Jf T = C™\t _1 (0), where r G C[x], is an open 
(quasi-)affinc variety of dimension n and the ring V(X T ) is identified with the 
localized Weyl algebra A n [r _1 ] := C^t" 1 ] <8>c[x] Ai(C). 

2.2. 25-modules. In general, the ring T>(X) carries a natural filtration 0(X) = 
Vq{X) C V\{X) G 2?2(X) G . . . provided by orders of differential operators (so 
that T>k is the subspace of operators of order < k) with [T>k,T>i] G T>k+i-i and 
T>(X) = Ufc>o£>fc(A"). The associated graded ring gr(V(X)) := ®k>o'Dk/'Dk-i is 
a commutative noetherian ring isomorphic to the affine algebra of the cotangent 
bundle of X, i.e. gr(V{X)) ^ 0(T*X). There is a natural mapping a : V(X) -> 
gr(2?(X)) which associates to a differential operator Q G T>(X) its principal symbol 
cr(Q), a polynomial function on T*X homogeneous along the fibers of T*X. The 
following properties are immediate: (a) cr(QiQ2) = o-(Qi)a(Q 2 ); (b) if 21 G £>(X) is 
a left ideal in 1>{X) then <r(2l) is an ideal in gr(T>(X)); moreover, (c) if Q\, . . . ,Qi G 
21 can be chosen in such a way that <r(Qi), . . . , cr(Qi) generate cr(2t) in gr(D(X)) 
then Qi, ... , Qi generate 21 in T>(X). As a consequence, we note that T>{x) is a left 
(right) noetherian ring. 

Let M be a left unitary D(X)-module. A V-filtration on M is a sequence M C 
Mi C . . . of C(A)-submodules such that (i) U fc > M fe = M and (ii) V k Mi C M fe+i 
for all G Z+. If the last inclusion (ii) holds with an equality for all sufficiently 
large I G Z + then the P-filtration is said to be good. In fact, M can be equipped 
with a good ©-filtration if and only if it is finite over V{X). Given a system of 
generators {ei,... ,e s } of a finite X>-module, one can define a good I?- filtration 
by setting := T>k{e\,... , e s ) for every k G Z + . Such Z?-filtrations are called 
standard, any two of them, say {M^} and {M' k }, being equivalent in the sense that 
M k C M^ +fco and Mj(. C M fe+fco for some k G Z+ and for all k G Z+. 

The important geometric invariant of a P-module M is its characteristic variety 
Vm- Let, for simplicity, X = C", and let M be a finite module on X over the Weyl 
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algebra V = A n (C). Consider the graded gr(2?)-module grp(M) := (Bk>oMk/Mk-i 
associated with some good Infiltration on M and take its left annihilating ideal in 
gr(X>): 

(2) Id (M) := {a £ gv(V) : a ■ gr v (M) = 0}. 

The variety Vm is then defined as the zero set of It>(M) in T*X = C 2n , i.e. 

(3) V M := {(x, 6 C 2n : a(x, £) = for all a £ I V (M)}. 

Clearly, Vm is a £-conic algebraic variety in C 2n whose "vanishing" ideal Jm is equal 
to the radical of I-p(M) (Hilbert's Nullstellensatz) and called the characteristic ideal 
of M. It turns out that both Vm and Jm depend only on M and not on the filtration 
chosen. 

If 7r : T*X — > X is the canonical projection, one defines the support of a T>- 
module M as a 7r-image of its characteristic variety 

(4) supp (M) := tt(V m ) = {x £ C" : 3(x, £) £ V M }. 

The dimension of Vm is called the Bernstein dimension d(M) of M. It is a funda- 
mental fact in the theory (cf. (l0[ that d(M) > n for every nonzero M. This 
suggests that one should single out the class of so-called holonomic 2?-modulcs 
B n (the Bernstein class) for which d(M) = n. Holonomic £>-modules possess a 
number of finiteness properties that make them of particular interest from the an- 
alytic point of view (see, e.g., p^|). For example, every M £ B n has finite length 
(i.e., any ascending chain of X>-submodules in M becomes stationary). In view of 
Stafford's theorem, this implies the cyclicity of M since T> is a simple ring of infinite 
length. As a whole, the Bernstein class B n is stable under certain basic operations 
on 2?-modules such as taking direct-inverse images, localization and the Fourier 
transform. 

In general, the precise determination of a characteristic variety is a difficult 
problem, even in the case of holonomic £>-modules. On the other hand, if a T>- 
module M has a characteristic variety Vm of a particularly simple structure one 
can expect to describe it explicitly. In Section |§|, we illustrate this remark in 
the simplest case: for holonomic £>-modules M £ B\ supported at a finite set of 
points in the complex plane, supp (M) — {Ai, . . . , \ s }. Using the basic operations 
on 2?-modules we then recover the Wilson construction of rank one commutative 
algebras of ordinary differential operators. The very existence of such algebras 
reflects a certain symmetry (or, to be more precise, "simplicity" ) of the underlying 
2?-module and its characteristic variety. 

3. Holonomic 2?-modules and Commutative Algebras of Ordinary 

Differential Operators 

In this section we review a construction of rank one commutative rings of dif- 
ferential operators in dimension n = 1 due to G. Wilson |37j. Our purpose is to 
analyze Wilson's scheme in terms of holonomic £>-modules over the Weyl algebra 
and its localizations. Reinterpreting Wilson's construction in these terms motivates 
the higher dimensional generalization to be presented in Section ^. 

3.1. 2?-modules with point support. Let X — C" and let M be a finite left T>- 
module on X with V = A n (C). Suppose that supp (M) = {A}, A £ C". In this case 
one can recover the structure of M explicitly. Indeed, it is clear that dime Vm < n, 
where Vm is the characteristic variety of M. In view of Bernstein's inequality this 
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implies d(M) = n and Vm — TTC n . Hence, M is a holonomic D-module, M G B n . 
If Mo C Mi C A/2 C . . . is a good nitration on M, the associated graded module 
gr-pM is supported at A. By the Nullstellensatz one has 

(5) mjMa)=0 

with some N G Z+ for any a(a) G grp(M), where ttia = C[x (x — X) is the (maximal) 
ideal in C[x] corresponding to the point A. But then, (51) is also true for any 
element in M: m^a = 0. To clarify the structure of M we start with the simplest 
example: A(A) := V/Vm x . Clearly, suppA(A) = {A} since m^ +1 Q E A(A) for 
every Q £ V of order ordQ < N. Denote by 5\{x) the canonical generator of A(A): 
5\{x) := 1 + T>m\ (so 8\ (x) is the image of 1 G V under the natural surjection 
T> — ► A(A)). It is easy to verify that 

A(A) = V(5\{x)) = C[d](6\(x)}, 

where C[d] is the ring of differential operators with constant coefficients in C n . 

The following lemma determines the algebraic structure of a finite £>-module 
with point support. 

Lemma 3.1. If M is a finite-type left T>-module on X = C" with suppM = {A}, 
A G C™, then M is isomorphic to a finite direct sum of copies of A(A) = V/Vm\. 
In particular, A(A) is a simple V-module (i.e. has no nontrivial submodules over 
V). 

For a proof of this result the reader is referred to |Q, Ch. 5 (see also [^9[ p0[). 

3.2. Left Ideals of the Weyl Algebra. Now let n = 1 and let M be any finite 
left P- module over the Weyl algebra V = Ai(C). We choose a system of generators 
{ei, . . . , e s }, M — P(ei, . . . , e s ), and fix the standard filtration Mi C Mi C M 2 C 
. . . on M associated with this system Mj~ := T>k(ei, . . . , e s ). Then, we denote by 
21 = 2l(ei, . . . , e s ) C T> the left ideal in T> that annihilates the elements {ei, . . . , e s } 

(6) 21 := Ann^ei,... ,e s ) = {Q(x,d) €V : Q(x, d) ej = 0, j = 1, . . . , *}. 
Put 21a_, := T>k n 21 and consider the associated graded ideal in gr(X>): 

gM2l) = 02l fc /2l fe _i 

fc>0 

containing the principal symbols of the differential operators in 21. Clearly, grjj(2l) 
is a homogeneous ideal in the polynomial ring C[a;,£] = gr(2?), and therefore it can 
be presented in the form 

(7) grp(2l) =©J fc (SQ£*. 

fc>0 

Here, Io(2t) C /1 (21) c /a (21) ■ ■ ■ is an increasing sequence of ideals in <C[x] that 
corresponds to leading terms of the differential operators in 21. Since C[x] is a 
principal ideal domain (recall, n = 1), each ifc(2l) is cyclic while the sequence {Ik} 
stabilizes: I q = I q +\ = Iq+2 = ■ ■ ■ starting from some q > 0. We write 7% = Tk(x) 
for the generators of Ife(2l) in C[x], and notice that each is determined uniquely 
(up to a constant factor). Let 

p := mm{fc G Z+ : 4(21) ^ 0} 

and 

5 := min {k G Z + : Zfe = // VZ > £;} 
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so that q > p. Then, the polynomials {t p (x)^ p , . . . ,T q (x)£ q } generate the graded 
ideal grp(2l) in C[x, £]. We refer to the integers p and a q := deg r q (x) as the order 
and degree of the ideal 21 respectively, and write Ord2l := p, Deg 21 := a q . A set 
of differential operators {L p (x,d), . . . ,L q (x,d)} C 21 such that <r(L r ) — r r (x)^ r , 
r = p, p+l, ... , q, is called a standard basis of the ideal 21 (cf. juj). The differential 
operator L p of the minimal order in a standard basis is determined uniquely, once 
a normalization of its leading coefficient is fixed. If follows from our previous 
discussion that a standard basis is indeed a system of generators of the ideal 21. 
In fact, it is not difficult to prove that 21 = T>(L p ) + T>(L q ) in agreement with a 
well-known result of Stafford that every ideal in the Weyl algebra A n (C) can be 
generated by two elements (see Q , Ch. I) . 

3.3. Wilson's Construction. Let M = ©^ =1 A(A J ) be a direct sum of finite left 
2?-modules with point support. It is easy to see that Vm = UT^.C and, hence, that 
supp M = {Ai, . . . , An}- We fix a C-linear finite-dimensional vector space V C M, 
homogeneous with respect to the direct decomposition above: V = ©jA(Aj) (~l V, 



and such that M = f (V). Then, according to Lemma 3.1, V has a finite basis of 
the form ej = pj (d) + Vm\j . Following Wilson |37| , we also consider the subring 
TZv C C[x] of polynomials P G C[x] with the property P.V C V, i.e. P. G Endc(V). 
Since supp M is finite, IZy is not empty. In particular, it contains all P G . 
for a sufficiently large k G Z + . The left ideal 21 := Annx>(V) annihilating V has the 
structure of a right module over TZ V . Indeed, (2l.P)V = 2l(F.V) C 21V = so that 
21 D 21.P by definition of 21. Note also that Ord2t = 0. 

Now we apply the operation of a Fourier transform M — > M v to the £>-module 
M. By definition, this is a (unique) bijection induced by the formal involution 
on the Weyl algebra T> = Ai(C): x \— * —d x , d x i— > x. One can prove |l^| that 
M v e Bi provided M e B\. If V v C M v is the Fourier image of V C M, then V v 
is annihilated by 2l v , the image of 21 under the Fourier involution in T>. The fact 
that Ord2l = implies that 2l v n C[d] is not empty. In particular, Deg2l v = 0. 

Let {G p (x,d), . . . , G q (x,d)} be a standard basis of the ideal 2l v so that 2l v = 
T>{G P ) +V(G g ). Clearly, we may choose G q € C[d], while G p (x, d) £ C[d] in view of 
the homogeneity of the space V. Recall that the operator G p is determined uniquely 
by the ideal 2l v (up to an inessential constant factor). We call the leading term of 
this operator r = t p (x) the t -function associated with the ideal 2l v . The idea now 
is to localize the Weyl algebra A\(C) at r = t p {x) and to consider the extension of 
the ideal 2l v in this localization. 

Thus, we set 

V T := Vlr- 1 } S Clx,^ 1 } ® c[x] A ± (C) 

and denote by 21^ the left ideal in T> T generated by 2l v . It follows that 21^ is a 
(principal) ideal with a unique generator D = D(x,d) := r(x)^ 1 o G p (x,d) G T> T , 
since ker G p (x,d) C keiG q (d) by construction. On the other hand, 21^ has the 
structure of a right module over the (commutative) ring IZy C C[d], the Fourier 
transform of IZv C C[x]. Hence, we conclude that for every L e TZy there exists 
an operator L G V T such that D o Lq = L o D in 21^ . The family of all differential 
operators L G T> T obtained in this way forms a commutative ring in the algebra 
V{X) of regular differential operators on X = C 1 \r~ 1 (0), isomorphic to TZ\>. It 
is not difficult to show (cf., |57j) that Speedy is a rational curve with cuspidal 
singularities located at the points of suppM. In this way, we recover Wilson's 
scheme of rational Darboux transformations in dimension n — 1 (cf . Q |2^] ) . 
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4. Algebraic Darboux Transformations for Partial Differential 

Operators 

The purpose of this section is to generalize Wilson's construction to higher di- 
mensions. We propose a procedure for generating (super)complete commutative 
rings of partial differential operators in arbitrary dimension whose common (for- 
mal) eigenfunction ip = iji{x,^) (being necessarily unique) has a particularly simple 
structure. By analogy with the one-dimensional case, we call this procedure the 
algebraic Darboux transformation. The fact that, in general, ij)(x,t;) has no trivial 
symmetry between spatial and spectral variables leads to a new example of the 
bispectral duality (in the sense of |l8|, |l9|, [57| ) in the case of several variables. More 
precisely, we extend the algebraic Darboux transformation to its "bispectral" coun- 
terpart and construct a commutative ring of (partial) differential operators in the 
spectral parameter sharing the same eigenfunction if) — iplx,^). However, unlike 
the one dimensional case, we are far from claiming that our procedure completely 
settles the bispectral problem in many dimensions. Our intention is to explore some 
hrst instructive examples in this direction. (For some other examples, see ||.) 

4.1. Commuting differential operators associated with £>-modules sup- 
ported by an algebraic hypersurface. 

4.1.1. The multiplicity free case. Let q{x) £ C[x] be an irreducible polynomial in 
n > 1 variables, x := (x\, . . . ,x n ), and let m q stand for the corresponding prime 
ideal in the polynomial ring. Extending the procedure developed in Section |[ we 
start with the 2?-module A(g) := T>/T>m q on C" supported at the zero set of q. 
We single out a cyclic submodule M := V(e) of A(q) generated by the element 
e := t(— &) + T>m q where r = t{x) is an arbitrary fixed polynomial in C[a;], and 
write 21 := Annp(e) for the left annihilator of e in T>. Clearly, 21 is a left ideal in 
the Weyl algebra of zero order, Ord2l = 0, since e = in A(q) for all sufficiently 
large N e Z+ (N > degr + 1). 

Let lZ e be the subring of differential operators in T> such that 

He := {Q G V : e G ker M (Q - A) for some A e C}. 

Put lZ e :— 1Z e nC[i]. Then, lZ e is not empty since 21 C 7Z e while 2lflC[a;] contains 
sufficiently large powers of the ideal m q . 

As in dimension n = 1, we apply to M the formal Fourier transform: M — * M v 
and similarly denote by 2l v , m^, TZ^ , . . . the objects dual to 21, m q , lZ e , . . . under 
the corresponding Fourier involution^]. Again, since Ord2l = 0, it is clear that 
2l v [~lC[<9] is not empty. Furthermore, 2l v is a right module over the (commutative) 
ring ■ 

Let X T be the open (quasi)-affine algebraic variety in C" obtained by removing 
the zero set of the polynomial r(x), i.e. X T := C"\r _1 (0). Then (see Exam- 
ple |, Sect. |2^), we can identify the ring T)(X T ) of regular differential operators 



on X T with the localization of the Weyl algebra at r = t(x), i.e., V T := T){X T ) 
C^T" 1 ] ®c[x] A n (C). Since T> C T> T , one can define the extension ideal 21^ of 
2l v in this localization, namely 2t^ := 2? r (2l v ). The following lemma clarifies the 
algebraic structure of 21^ . 

Lemma 4.1. 21^ is a principal left ideal in T> T . 



1 Here, m,j, 7?. e ,. . . are regarded as being imbedded in the Weyl algebra A n (C). 
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Proof. By construction, 21 is a left ideal in T> annihilating the principal generator 
e v := t{x) + Vm y q of the 2?-module M v . Fix N G Z + such that iV > degg, and 
define G(x,d) := t n (x) o o t(x)- 1 G X> T . Then, it follows that G(x,d) G 2l v . 
Indeed, in view of the Weyl algebra commutation relations one has ad^ [q(d)] = 0, 
where ad^ stands for the iV-iterated commutator with r = t(x). Since 

ad? (x) [q(d)} = £(-l) fc ( W)""* o q(d) o ,-(*)*, 
fc=0 ^ ' 

the expression t(x) o q(d) is divisible by t(x) from the right within A n (C). On 
the other hand, it is obvious that G(x, d) o t(x) = (mod Pra^). 

Put D = D(x, d) := t(x)' n+1 oG(x, d), so that D G 21^ in view of the definition 
of 21^. Then, we claim that the ideal 21^ is generated by D in V T . Indeed, if 
Q(x,d) G 21^, then Q is factorizable as Q(x,d) = S(x,d) oT(x,d) with T G 2l v 
and S G V T . By definition of 2l v , we have T(x,d) o r(x) = (mod £>m^), and 
hence T(x, d) o t(x) — T(x, d) o q(d) for some T(x, d) G V. Multiplying both 
sides of the last equation by S(x, d) on the left and by t(x)~ 1 on the right, we get 
Q(x, d) = f'(x, d) o D(x, d) where f' = fo t(x)' 1 G V t . Thus, 21^ is a cyclic left 
ideal in V T . □ 

Now we are ready to state the main theorem which completes our construction 
of algebraic Darboux transformations related to the P-module A(q). 

Theorem 4.2. Associated with M C A(q) there exists a commutative ring TZm C 
T> T of regular differential operators on X T = C n \r (0) isomorphic to TZ e . If K(q) 
and A(t) are the complex linealitiesf\ of the polynomial q and r and if A(q) U A(r) ^ 
C™ , the ring TZm is nontrivial in the sense that TZm cannot be transformed into 
7Z y C C[d] by conjugation by a regular function on X T or through a local nonsingular 
change of variables. 

Proof. By definition, 21^ is generated by 21 v as a left module over V T . Hence, it 
inherits (from 2t v ) the structure of a right module over the commutative ring TZ y . 



Therefore, in view of Lemma 4.1, we conclude that for every Lq G TZ y there exists 
an operator L G T> T such that D o Lq = L o D in 21^ with D = t(x) o q(d) o t(x)^ 1 . 
Denote by TZm the set of all differential operators L G T> T constructed through this 
procedure. Then, TZm forms a commutative subalgebra in the ring T>(X T ) since the 
latter has no zero divisors. When A(q) U A(t) ^ C™, A ± (g) n A (r) is not empty 
and [q(d), t(x)) 0, so that D C[d]. This implies that TZm is neither equal nor 
conjugate to TZ y e . Moreover, since the operators in TZm have constant principal 
symbols, they cannot be reduced to their principal parts through a local change of 
coordinates. □ 



Remark 4.3. Recall that, by construction, TZ e (and hence TZ^) is not empty. In 
particular, TZ e contains the powers D m^ +1 D ... of the ideal m q starting 
with N = degr(x) + 1. By Fourier duality, this implies that (m^)^ C TZ^ for all 
N > degr(x) + 1. In fact, if Lq G (m^ ) N , then the corresponding operator L G TZm 



2 By a complex lineality of a polynomial P S C[ir] we mean the maximal linear subspace A(P) 
of C n such that P can be restricted to a polynomial on the quotient C n /A(P). 
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can be constructed explicitly. Indeed, we have Lq = p(d) o q(d) N with p(d) G C[d]. 
On the other hand, in view of the Weyl algebra commutation relations, the identity 

«< 8) [r(ar)] := [q(B), [q(d), . . . [q(d),r(x)} ...]]= 0, 

or more explicitly, 



jr(-i) k ( N k )<i(d) N - k oT( x )o q (d) k 

k=Q ^ ' 



= 



holds for all N > degr + 1. Whence, we get 

N / \ 

q (d) N o r(x) = $>l) fc+1 (Z) l(df- k o T o q(df, 

k=l ^ ' 

so that Lo can be factorized in the form 

(8) L = Q{x,d)oD(x,d). 
Here 

N 

(9) Q(x, d) := ^(-l) fc+1 ( * )p(d)q(d) N - k o r(.x) o g(S)*- 1 o r^)" 1 

fe=i ^ ' ' 

and 9) := t(x) o g(<9) o t(x) _1 . By interchanging the factors in (||), we obtain 
the corresponding operator in TZm'- 

(10) L := D(x,d) o Q(x,d). 

This procedure is a precise analog of the classical one-dimensional Darboux trans- 
formation p7| . 

4.1.2. TVie general case. This procedure of algebraic Darboux transformations can 
be extended to the case of arbitrary multiplicities. Again, we start with a pair 
of complex polynomials q — q(x) and r = t(x) in n > 1 independent vari- 
ables (without assuming q to be primitive). In this case, we can decompose 
q = Cq^ 1 q% 2 . . . q^ s into a product of powers of pairwise distinct irreducible poly- 
nomials qi = qi{x) G C[x] with (strictly) positive integer multiplicities k t £ Z + , 
C G C x , and we write := m Qi for the corresponding prime ideals in C[x]. To 
each ki we then associate the C-linear vector space Vi C C[x] spanned by all mul- 
tiples of t(x) by polynomials of order less than ki, i.e. we set 

Vi = spanc(a(x)r(x) G C[x] : dega < ki). 

With the notation of Section [4.1.l| we define A(gi, ... ,q s ) := A(gj) as a left 

X>-module over the Weyl algebra on C n and consider its submodule M generated 
by the elements a :— pi(—d) + Vrxii with all pi 6 Vi, i = 1, . . . , s. 

Let M v be a Fourier transform of the D-module M. Clearly, M v = V(e\ , e%) , 
where e\ := Vi + L>m^, while stands for the image of mi under the (formal) 
Fourier involution in V. We write 2l v C T> for the left ideal in T> annihilating 
the elements e^ , i = 1, . . . , s, 2l v := Annp(e^ , . . . , e^f ) and consider its extension 
21^ := 2? r (2l v ) in the localization of the Weyl algebra V T at r. 



In this setting, we can refine the argument leading to the proof of Lemma 4.1 



i.e. we can show the cyclicity of the ideal 21^ within V T . Indeed, if Q(x,d) G 21^, 
one has Q(x, d) = S(x, d) o T(x, d) with S(x, d) G V T and T(x, d) G 2l v C V. 
The latter implies that T(x,d) o p(x) — T p (x,d) o for any p(x) G Vi and 
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some T p £ T>. Since p{x) = a{x)r{x) with a(x) £ C[x], dega < fc,, we can 
rewrite the above equation in the form R(x,d) o a(x) = T p (x,d) o qi(d), where 
R(x, d) := T(x, d) o t(x) £ T> T . It follows that all the derivatives of the full symbol 
R(x, £) of the operator R up to order hi — 1 with respect to £ vanish on the zero set 
of the ideal in C". Elementary induction in fcj with the use of the Nullstellensatz 
gives immediately that R(x,£) = Ri{x, ^)q i i (£) and, hence, T{x,d) is divisible on 
the right by qi(d) ki ot(x)^ 1 for all £ = 1, 2, . . . , s. Finally, due to the commutativity 
of the operators Di := t(x) o qi(d) ki o t{x)~ 1 £ V r , we conclude that T must be of 
the form T(x, d) = T'(x, d) o D(x, d) where D := Di = t(x) o q(d) o r(x) _1 and 
T" £ P T . In this way, the operator 9) is a principal generator of the ideal 21^ 
in 2? T . 

Let 7^ e := {P G C[x] : P . spanc(ei, . . . , e s ) C spanc(ei, . . . e s )} and let TZ^ be 
its formal Fourier transform. Both TZ e and 7£g are not empty: they contain in par- 
ticular the ideals flm^* and n(m^) fci respectively. Then, the left ideals 2l v and 21^ 
have the structure of right modules over the commutative ring TV e , and therefore 



we again arrive at the result stated in Theorem 4.2. Namely, we obtain a commuta- 
tive ring TZm C V t of regular differential operators on X T = C"\r _1 (0) isomorphic 
to iZg . This completes the procedure of algebraic Darboux transformations in the 
case of arbitrary multiplicities. 

4.2. The Bispectral Duality. The formal spectral problem for the commuta- 
tive ring TZm has a unique solution of a particularly simple structure. Indeed, 
by construction, the (quasi-)regular function on T*X T given by if) = %l>(x,£) := 
D{x,d)e^ x '^ is a common eigenfunction for all differential operators in TZm- We 
call ip the Baker- Akhiezer function associated with TZm- 

The following elementary argument shows that there exists a certain complete 
commutative ring TZ M of partial differential operators on X q = C™\g _1 (0) written 
in terms of the spectral variable £ for which ij) is also a (unique) common eigen- 
function. 

Put N = &egq(x) and denote by P the principal ideal in the polynomial ring 
C[x] generated by t(x) n . For every 7 £ I b we have r f{x)iji{x, £) = Qy(x,d x )e^' 
with some Q 1 £ A n (C). Then, there obviously exists another differential opera- 
tor Qj($,,d^) in A n (C) written formally in terms of £ such that Qj(x, d x )e( x & = 
Q\(£,d 6 )e (x ^\ and hence 

(11) Cb(ZA)e M =7(#,5). 
On the other hand, we have 

(12) r(dt)1>(x, = r(x)q(d x )e^ = q^) T (x)e^ 



in view of the definition of ip. It follows immediately from (|ll| ) and ( |12[ ) that 
(13) L 7 (& ds)tl>(x, = 7 (x)r(a#(z, 0. 

where d%) := Qt(£, d%) o g(^) -1 o r(d^). The family of all operators dg) 

constructed in this way forms a commutative ring TZ M of partial differential oper- 
ators in the localized Weyl algebra V q = C[£, 9(C) 1 ] ®C[£] -^n(C) isomorphic to I b . 
This ring can be regarded as a bispectral dual of TZm in the sense of [ |l8| p7| . The 
simple procedure we have applied to generate the ring TZ M is, in fact, an extension 



of the bispectral Darboux transformation scheme in n = 1 as discussed in 22 
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In conclusion, we note that the ring TZ\ I: though complete, is not necessarily 
a maximal commutative subring in T> q , i.e., the centralizer Z(7Z. M ) of TZ^ M may 
happen to be larger than 1Z M itself. It is an interesting problem to clarify the 
structure of Z{Ti}' M ) in general. 



5. Examples, Applications and Concluding Remarks 

The purpose of this section is to display some explicit examples of the theory 
presented so far and to illuminate some of its implications. 

Let n = 2. We choose r(xi,x 2 ) = x\ — x 2 and g(xi,x 2 ) — xix 2 — A, A e C. If 



A ^ 0, the polynomial q is irreducible. According to Theorem 4.2, for any constant 
coefficient operator Q 6 C[di, <9 2 ] and Lq — (8182 — A) 3 there exists a unique partial 
differential operator Lq satisfying the identity 

(14) DoQoL = L Q oD, 

where 



(15) D:=T(xi,x 2 )°q{di,d 2 )oT(xi,X2) 1 



di-^A (d 2 H — j~ ) ' x - 

xf - X 2 J \ X\-X2, 



For example, if Q = 1, formulas (p|)-(pij|) (see Remark p~3|) give explicitly 



{xf-x 2 ) \X\-X2) 

-3 (-4a;i + \x\ - 2Xx\x 2 + Xxf) „ „ -6 (13a;? + 5x 2 ) n n „ 

-dfdi + ^- ^-didi 



(x\ ~ x 2 ) ~ (x\ - x 2 y 

6 (18a;? + Ax? + 30xix 2 - Axf x 2 - \x\xl + Xxf) „ -3 „, „ 

— 1 -82 H j "1^2 

(x 2 -x 2 ) {x\-x 2 ) 

54xi o2 „ 3 (6A 2 x?x 2 - 4A 2 x 2 x 3 + A 2 x|) o „ 
+ -5 ^dfch + -± y ^8,82 

(x\-X 2 ) (x\-X 2 ) 

3 (-118a; 2 - 4Axf + A 2 x? - 26x 2 + 8Ax?x 2 - 4A 2 xfx 2 - AXx^fj „ „ 

8\8 2 



{x\ - x 2 ) 

6 (I36xf + 9Axf + 104xix 2 - 15Axfx 2 + 3Ax 2 x| + 3Aa;f ) -3 3 

n — 7 " 2 + 71 — 

(xf-x 2 ) (x(-x 2 ) 

3 (24xi + Xxf - 2Ax?x 2 + Ax|) 2 

^1 



(x 2 - x 2 ) 

-12 (52x 2 + 3Axf + 8x 2 - QXx\x 2 + 3Axix|) 



di 



(x\ - x 2 ) 

1920x? - A 3 x} 2 + 960xix 2 + 6A 3 x}°x 2 - 15A 3 xf x 2 + 24Ax 3 
(xj - x 2 ) 

-A 3 x6 + 20Ax^ (6 + A 2 x 3 ) + %Xx\x\ (12 + X 2 xf) - 3Axtx 2 (72 + 5A 2 x 3 ) 

(xf - x 2 ) 
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Note that the operators Lq are nontrivial in the sense that they cannot be reduced 
to the constant operators Lq := Q o Lq by means of conjugation or change of 
variables. 

The simplest operator in the spectral parameter which has the same eigenfunc- 
tion as the Lq's can also be computed in an explicit form using (pi]), ( O ) and (|l3|). 
Namely, we have 

(17) L\^d i ) = Q b 1 oq^ 1 ^ 2 )- 1 or(d i ), 

where 

Qifo % = (66 - A)^ - 26^ - 2(66 - A)££% 

+ Zidl + (66 - X)dl + 2&%% - 6% - 4% . 

The next example reveals the link to the theory of Huy gens' principle on Minkowski 
spaces. Let n be an arbitrary even number, n > 4. Take t(x) — x\ — x 2 and 
q(x) = x\x 2 — x\ — x\ — ■ ■ • — x^- Clearly, g is again an irreducible polynomial and 



we can apply Theorem 4.2. The simplest operator in the commutative ring TZm has 
the same form as (|l6|), where the parameter A is formally replaced by the Laplacian 
A„_2 = 9| + 8\ + ■ ■ ■ + d„ in n — 2 extra variables. This makes sense since A enters 
polynomially in to Lq and L\, The operator Lq = (8182 — 8\ — ■ ■ ■ (9^) 3 is a cube 
of the classical wave operator, and hence, it satisfies Huygens' principle for all even 
n > 8 (cf., e.g., H). 

Adapting the technique developed in H, we can prove that the operator L\ 



d\d\-\-. . . obtained by the formal substitution A A„_2 in (16) is also a nontrivial 
Huygensian hyperbolic operator in even dimensions n > 12. Note that L\ cannot 
be presented as a power of the second order wave-type operator on the Minkowski 
space M™. Therefore, we get a new interesting example related to the classical 
Hadamard's problem (see [pL EG] and references therein). 
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